Weak ferromagnetism and field-induced spin reorientation in K2V3O8 
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Magnetization and neutron diftraction measurements indicate long-range antiferromagnetic ordering 
below Tm=^ K in the 2D, S=l/2 Heisenberg antiferromagnet K2V30g. The ordered state exhibits 
"weak ferromagnetism" and novel, field-induced spin reorientations. These experimental obser- 
vations are well described by a classical, two-spin Heisenberg model incorporating Dzyaloshinskii- 
Moriya interactions and an additional c-axis anisotropy. This additional anisotropy can be accounted 
for by inclusion of the symmetric anisotropy term recently described by Kaplan, Shekhtman, Entin- 
Wohlman, and Aharony. This suggests that K2V3O8 may be a very unique system where the 
qualitative behavior relies on the presence of this symmetric anisotropy. 
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One of the most exotic and intriguing anisotropies 
which can occur in magnetic solid state systems is the 
antisymmetric Dzyaloshinskii-Moriya (DM) interaction. 
These interactions, originally proposed by Dzyaloshin- 
skii 0] and later derived by Moriya I through extension 
of Anderson's superexchange model ||], arise from spin- 
orbit coupling in the presence of sufhciently low local lat- 
tice symmetry. The form of this interaction is typically 
written D • (Si x S2) where D is a vector determined 
by the symmetry of the lattice. Competition between 
the vector product in the DM interaction and the scaler 
product in the Heisenberg interaction often leads to novel 
phenomena such as spin-canting [Q and spiral spin struc- 
tures ||-|. 

Relatively recently, Kaplan and later Shekhtman, 
Entin-Wohlman, and Aharony (KSEA) (moti- 
vated by relevance to parent compounds of high-Tc su- 
perconductors), re-examined the DM interaction with 
particular emphasis on an often ignored symmetric com- 
ponent. This symmetric anisotropy always accompanies 
the typical antisymmetric DM term but is typically ig- 
nored as its strength is an order of magnitude weaker. 
KSEA showed that, in the general case where there is 
no frustration, this symmetric term cannot be ignored 
as its inclusion restores the original 0(3) symmetry of 
the Heisenberg Hamiltonian. We report in this letter 
the novel properties of a S=l/2, 2d Heisenberg system, 
K2V3O8, and suggest that they can be explained by a 
model including both symmetric and antisymmetric DM 
interactions. This implies that K2V3O8 is a novel system 
where the qualitative behavior relies heavily on inclusion 
of the symmetric component of the DM interaction. 

K2V3O8 crystallizes in a tetragonal unit cell with 
space group P4bm and lattice constants a=8.870 A and 
c=5.215 A The structure is shown in Fig. 1(a) and 
consists of magnetic V^^'^-Os pyramids and non-magnetic 
V^"'"-04 tetrahedra with interstitial K+ ions. Previous 
magnetic measurements consisted of powder magnetiza- 
tion from 5-300 K [Q which were best described by a 2D 
Heisenberg model with coupling constant J=12.6 K and 



a g- value of 1.89. EPR measurements on single crystals 
suggest very small anisotropy with g- values of gc= 1.922 
and gaf,= 1.972 ||l^]. Liu ct al. suggested from ob- 
servation of field-dependent magnetization that K2V3O8 
may order at lower temperatures. 

Single crystal plates of K2V3O8 (typical dimensions: 
1 X 1 xO.l cm^) were grown in a platinum crucible by cool- 
ing appropriate amounts of VO2 in a molten KVO3 flux. 
Additional details of the crystal growth will be reported 
elsewhere. ||l^ One of these single crystals was mounted 
in a Quantum Design SQUID magnetometer equipped 
with a sample rotator and the magnetization is shown 
in Fig. 2 as a function of both temperature (a) and ap- 
plied magnetic field (b) for fields applied along both the 
c-axis and within the tetragonal basal plane. The M(T) 
data (Fig. 2(a)) were measured by cooling in the pres- 
ence of an applied magnetic field while the M(H) data 
(Fig. 2(b)) were taken by cooling to base temperature in 
zero-field. From Fig. 2(a), the magnetization is seen to 
be very isotropic for temperatures in excess of about 8 
K confirming the Heisenberg nature of the interactions. 
At high temperatures, we observe a Curie- Weiss suscep- 
tibility with an effective moment (/i=g/iB[S(S-)-l)]^/^) of 
1.7/iB per formula unit (FU). This result is consistent 
with the presence of two non-magnetic V-ions per FU 
and one magnetic ion with S=l/2. 

The temperature dependence in the presence of a weak 
magnetic field shows a clear ordering phase transition be- 
low a temperature of about 4 K. The rapid decrease of 
the magnetization below this temperature with field ap- 
plied along the c-axis is indicative of antiferromagnetic 
ordering state with an easy axis along c. However, mea- 
surements within the basal plane show a clear ferromag- 
netic enhancement upon entering the ordered state. The 
ferromagnetic ordered moment is small (~10~^ /is) and 
thus we conclude that the system is antiferromagnetic 
with spins primarily along the c-axis and a small canted 
moment in the basal plane. This represents the first ob- 
servation of an ordering phase transition in this material. 

The behavior of the magnetization as a function of 
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FIG. 1. (a) The crystal structure of K2V3O8 composed of 
magnetic V*"''-05 pyramids and non-magnetic tetra- 
hedra. (b) The projection of the V*^ positions showing the 
location of the magnetic moments. The red arrows represent 
the zero-field spin configuration and the blue arrows denote 
the behavior of the system in the presence of a basal plane 
magnetic field (in the direction shown by H). 



field for temperatures well below the Neel temperature is 
shown in Fig. 2(b). For field applied along c, we observe 
an abrupt increase in magnetization for fields in excess 
of about 0.85 T consistent with a discontinuous spin-flop 
phase transition, further confirming the easy nature of 
the c-axis. For fields within the basal plane, we observe 
an intriguing field-induced phase transition which occurs 
at a magnetic field of H5ij=0.65 T. There was no evidence 
of irreversibility in M(H) for either field orientation or in 
angular rotations of the sample in the presence of an 
applied field and no evidence of thermal hysteresis was 
observed. 

To investigate the detailed microscopic spin arrange- 
ment, neutron diffraction measurements were performed 
on the HBlA and HBl triple axis spectrometers at the 
High Flux Isotope Reactor in Oak Ridge National Lab- 
oratory. For the measurements in zero-field, the sample 
was mounted in the (hOl) scattering plane and compari- 
son of scattering at 1.6 K and 10 K indicated the presence 
of magnetic Bragg peaks (in the (hOl) plane) described 
by integer indices with h odd. These peaks suggest a 
magnetic unit cell in which the corner and face-center 
moments within the basal plane are aligned antiparallel 
with parallel alignment along the c-axis. Examination 
of the integrated intensities of 6 magnetic reflections in- 
dicated that the moments lie along the c-axis consistent 
with conclusions drawn from the magnetization measure- 
ments. The value of the ordered moment (/i=g/i_BS) was 
found to be 0.72(4) at T=1.6 K, reduced from the 
expected value of 1 ^b- This can likely be attributed to 
quantum fluctuations in the low-dimensional S=l/2 sys- 
tem. The zero-field magnetic structure is shown schemat- 
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FIG. 2. (a) The temperature dependence of the DC sus- 
ceptibility in a small applied magnetic field. For tempera- 
tures in excess of about 8 K, the system is seen to be very 
isotropic. Below 4 K, the system undergoes a phase transition 
to an antiferromagnetic state characterized by an easy c-axis 
and weak ferromagnetic behavior in the basal plane, (b) The 
magnetization as a function of applied magnetic field for tem- 
peratures well below the ordering temperature. Field-induced 
phase transitions are observed for field applied in the basal 
plane (Hs_r) and along the c-axis (Hsf). 



ically by the red arrows in Fig. 1(b). 

To examine the behavior in the presence of a basal 
plane magnetic field, the crystal was mounted in the 
(/lOZ) scattering plane and a vertical field applied along 
the (010) direction. The field dependence of the (100) 
and (101) magnetic Bragg peaks are shown in Figs. 3(a) 
and 3(b). As the intensity of magnetic neutron scatter- 
ing is proportional to the component of magnetic mo- 
ment normal to the wavevector transfer, Q, the observed 
disappearance of the (100) magnetic Bragg peak and the 
concomitant enhancement of the (101) reflection are con- 
sistent with a continuous rotation of the spins from the 
(001) to the (100) direction with increasing magnetic field 
strength. The rotation is complete at a field strength of 
H5fl=0.65 T in excellent agreement with the anomaly in 
the magnetization. This is shown schematically by the 
blue arrows in Fig. 1(b). Thus, we conclude that this 
transition is a peculiar spin-reorientation phase transi- 
tion whereby the spins rotate from the c-axis to the basal 
plane while remaining normal to the applied field direc- 
tion. 

A second crystal was mounted in the (hkO) scatter- 
ing plane to allow measurements with field applied along 
the c-axis. The results of these measurements are shown 
in Fig 3(c) and (d). As a function of applied magnetic 
field, the (100) Bragg peak intensity drops abruptly upon 
passing through a field of 0.85 T suggesting a discontin- 
uous reorientation of the spins. The temperature depen- 
dence in an applied field of 1 T shows a strikingly dif- 
ferent behavior for the (100) and (010) Bragg peaks and 
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FIG. 3. The background subtracted intensity of diffracted 
neutrons as a function of applied magnetic field (along the 
(010) direction), for temperature well below the Neel temper- 
ature, is shown for the (100) and (101) magnetic Bragg peaks 
in panels (a) and (b). The solid lines represent calculations 
based on the model described in the text, (c) The field de- 
pendence of the (100) Bragg refiection with field appHed along 
the c-axis, for temperature well below the ordering tempera- 
ture, showing a discontinuous spin-fiop phase transition, (d) 
The temperature dependence of (100) and (010) in an apphed 
field of 1 T. 



the weaker intensity for (010) suggests that the moments 
have "flopped" from the c-axis into the basal plane se- 
lecting a specific direction within the basal plane in close 
proximity to the (010) direction. No evidence for a struc- 
tural distortion has been observed and the selection of a 
specific direction likely results from a slight misalignment 
of the magnetic field with the c-axis. This misalignment 
was less than 1° for the measurements performed but 
even this level may induce a net field in the basal plane 
of sufficient strength to overcome domain energies thus 
selecting a specific magnetic domain. Such a domain se- 
lection has been observed in the tetragonal antiferromag- 
net MnF2 0- 

As mentioned in the introduction, weak ferromagnets 
often have at their origin the competition between the an- 
tisymmetric DM interaction and the symmetric Heisen- 
berg interaction. To consider the DM interactions in 
K2V3O8, we will consider the projection of the crystal 
structure onto the tetragonal basal plane as shown in Fig. 
4(a). The plane of inversion symmetry between near- 
neighbor V'*"'' spins results in DM vectors located in that 
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FIG. 4. (a) The projection of the crystal structure onto 
the tetragonal basal plane. The magnetic moments are lo- 
cated at the center of the shaded squares and the triangles 
represent the V^'''-05 tetrahedra. The components of the 
Dzyaloshinskii-Moriya vectors {Dxy and Dz) are shown in the 
figure relative to the central spin, (b) The calculated compo- 
nent of magnetic moment along the field direction as a func- 
tion of applied field for field applied along both a and c. This 
result compares very well to the M(H) data plotted in Fig. 
2(b). 

inversion plane Q giving components along (110) and 
(001) directions, denoted by Dxy and Dz respectively 
(these directional vectors are indicted in the figure rela- 
tive to the central site) . The four- fold rotation symmetry 
through the central site gives rise to the configuration of 
D vectors shown in the figure. From the neutron diffrac- 
tion measurements, we know that the magnetic unit cell 
contains two spins in both zero and non-zero fields and, 
hence, we will consider a two-spin, classical Hamiltonian. 
Noting that in averaging over the the near-neighbor spins 
within a unit cell, the components of the DM vectors in 
the basal plane (i.e. the Dxy terms) cancel, we write the 
ground state energy as 



F = JMi • M2 + 13^ [Ml X Ma]^ 



(1) 



where Mi and M2 are the sublattice magnetizations 
and [Ml x M2]2 is the z-component of the vector cross- 
product. 

This ground state energy will be minimized when the 
spins are in the basal plane and, consequently, an ad- 
ditional c-axis anisotropy is required to account for the 
easy c-axis observed experimentally. Thus, we add a c- 
axis anisotropy of the form EzMizM2z and, with addition 
of the Zeeman energy for nonzero field, we arrive at the 
ground state energy. 



F = JMi ■ M2 + Dz[Mi X M2]. 
-I- 5/isH • [Ml +M2]. 



E.MizAh 



(2) 



This rather simple, classical ground state energy de- 
scribes all the low temperature properties of the system. 
For appropriate values of and , the spins will point 
along the c-axis with canting observed in small fields due 
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to the normal competition between DM and Heisenberg 
interactions. For fields applied in the basal plane, the 
Zeeman energy and the DM interaction favor arrange- 
ment of the spins in the basal plane. This competes with 
the c-axis anisotropy causing a continuous rotation of the 
spins to the basal plane with increasing field strength. 
The spins remain primarily normal to the applied field 
direction to minimize the Zeeman energy. Finally, with 
field applied along the c-axis, a spin-flop transition oc- 
curs due to competition between Zeeman energy and the 
anisotropics in the system causing an abrupt rotation of 
the spins into the basal plane. 

As we have two field-induced phase transitions, we can 
estimate the values for D^/J and E^/J required to pro- 
duce phase transitions at the measured field values (using 
J=12.6K and g=1.89 as determined from Liu et al. @|). 
This results in a value of D^/J=0.22 and E^/J=3.8x 10"^ 
(this value of D/J seems large considering the measured 
g-values HJlll but compares favorably to the DM spi- 
ral system BaCuGe207 which has a similar basal plane 
structure and for which D/J was found to be ^0.18 [^). 
Using these values, we can calculate the induced mag- 
netic moment along the field direction for field applied 
along both the a- and c-axes and the results are plotted in 
Fig. 4(d) as a function of applied magnetic field in Tesla. 
These results can be directly compared to the M(H) data 
plotted in Fig. 2(b) and the qualitative agreement is 
remarkable suggesting that this relatively simple model 
adequately describes the low temperature properties of 
the system. To further emphasize this, we have super- 
imposed on Fig. 3(a) and 3(b) the calculated magnetic 
Bragg peak intensity for the (100) and (101) reflections 
in the presence of an applied magnetic field along the 
(010) direction. These results are represented by the solid 
lines in Figs. 3(a) and 3(b) and the agreement between 
the calculations and the data is excellent, particularly for 
the (100) reflection. The agreement for the (101) reflec- 
tion is very good for field strengths up to H5 a but seems 
to deviate for higher fields. The reason for this disagree- 
ment is unclear but could result from a fleld-dependent 
<S>. Such a scenario can occur via field-dependent sup- 
pression of quantum fluctuations, as has been observed 
in some ABX3 antiferromagnets [p^ . 

To shed some light on the nature of the c-axis 
anisotropy, we calculate the symmetric anisotropy term 
of KSEA. The addition of this term causes a rescaling of J 
to (J-D'^/4J) and the inclusion of a symmetric anisotropy 
term (l/2J)Mi-A-M2 where A is a 3x3 matrix with el- 
ements A„„=D„D^ ]ri| , |l9t . Making these changes to the 
Hamiltonian and averaging over near-neighbor spins pro- 
duces the ground state energy, 

F=iJ-^)M^■M2 + D, [Ml x M2]. 

+ i^-^)Mi.M2.+gfiB'ti-[Mi+M2]. (3) 



This energy has precisely the same form as Eq. 2 and, 
consequently, the c-axis anisotropy falls out naturally 
upon inclusion of both the symmetric and antisymmet- 
ric anisotropics. The maximum value for the c-axis 
anisotropy in Eq. 3 is D^/2J which, for D^/J of 0.22, 
gives a maximum anisotropy (E^/J) of 2.4x10"^^ Within 
the level of approximation involved in this model, this is 
in reasonably good agreement with the above value of 
3.8x10^^. Consequently, we conclude that K2V3O8 rep- 
resents a unique system where the qualitative features 
rely heavily on the inclusion of the KSEA symmetric 
anisotropy term. The only known, clear evidence of this 
symmetric anisotropy occurs in the DM spiral system 
Ba2CuGe207 where the addition of this term is neces- 
sary to produce quantitative agreement between exper- 
iment and theory As a caveat, it is important to 
note that we cannot rule out other possible mechanisms 
for this additional anisotropy. Rather, we simply note 
that inclusion of this symmetric anisotropy alone, which 
is necessary to properly consider DM interactions, seems 
to adequately describe the properties of the system. 

In summary, we have observed an ordering phase tran- 
sition in K2V3O8 with a primarily antiferromagnetic or- 
dered state accompanied by weak ferromagnetism. Field- 
induced spin reorientation phase transitions have been 
observed with fleld applied along both the c-axis and 
in the tetragonal basal plane. The spin arrangement 
has been determined by single-crystal neutron diffrac- 
tion measurements. Remarkably, this rich magnetic be- 
havior can be well described by a simple, classical, two- 
spin Heisenberg model with inclusion of DM interac- 
tions and an additional c-axis anisotropy. This additional 
anisotropy can, at least partially, be accounted for by in- 
clusion of the KSEA symmetric anisotropy term making 
K2V3O8 a unique system where introduction of this inter- 
action is necessary to describe the qualitative behavior 
of the system. 
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